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We discuss the phase diagram of moderately dense, locally neutral three- 
flavor quark matter using the framework of an effective model of quan- 
tum chromodynamics with a local interaction. The phase diagrams in the 
plane of temperature and quark chemical potential as well as in the plane 
of temperature and lepton-number chemical potential are discussed. 
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2 Riister, Werth, Buballa, Shovkovy, and Rischke 

1. Introduction 

Theoretical studies suggest that baryon matter at sufficiently high density 
and sufficiently low temperature is a color superconductor. (For reviews 
on color superconductivity see, for example, Ref. n) In nature, the highest 
densities of matter are reached in central regions of compact stars. There, 
the density might be as large as lOpo where po ~ 0.15 fm~^ is the satura- 
tion density. It is possible that baryonic matter is deconfined under such 
conditions and, perhaps, it is color-superconducting. 

In compact stars, matter in the bulk is neutral with respect to the 
electric and color charges. Matter should also remain in (3 equilibrium. 
Taking these constraints consistently into account may have a strong effect 
on the competition between different phases of deconfined quark matter at 
large baryon densities.^ 3 4 5 6 

The first attempts to obtain the phase diagram of dense, locally neutral 
three-flavor quark matter as a function of the strange quark mass, the quark 
chemical potential, and the temperature were made in Refs. ElandEI This 
was done within the framework of a Nambu-Jona-Lasinio (NJL) model. It 
was shown that, at zero temperature and small values of the strange quark 
mass, the gr oun d state of matter corresponds to the color-flavor-locked 
(CFL) phase.EElAt some critical value of the strange quark mass, this is 
replaced by the gapless CFL (gCFL) phase. ^ In addition, several other 
phases were found at nonzero temperature. For instance, it was shown t hat 
there should exist a metallic CFL (mCFL) phase, a so-called uSC ph ase^^ 
as well as the standard two-flavor color-superconducting (2SC) phas^^^^l 
and the gapless 2SC (g2SC) phase.l^ 

In Ref. T, the effect of the strange quark mass was incorporated only 
approximately through a shift of the chemical potential of strange quarks, 
/is — > /is — mg/(2/i). While such an approach is certainly reliable at small 
values of the strange quark mass, it becomes uncontrollable with increasing 
the mass. The phase diagram of Ref. [7| was further developed in Refs. El 
and 1141 where the shift-approximation in dealing with the strange quark 
was not employed any more, although quark masses were still treated as 
free parameters, rather than dynamically generated quantities. In Refs. 1151 
and [16' the phase diagram of dense, locally neutral three-flavor quark mat- 
ter was further improved by treating dynamically generated quark masses 
self-consistently. Some results within this approach at zero and nonzero 
temperatures were also obtained in Refs. ISl IT71 [THI and [T^ 

The results of Refs. El and El Etre presented here. Only locally neu- 



February 4, 2008 19:55 



WSPC/Trim Size: 9in x 6in for Review Volume 



chapter'vl 



Phase diagram of neutral quark matter at moderate densities 3 

tral phases are considered. This excludes, for example, mixe and 
o^ 

crystalline phases. Taking them into account requires a special treat- 
ment. We also discuss the effect of a nonzero neutrino (or, more precis ely, 
lepton-number) chemical potential on the structure of the phase diagram.l^ 
This is expected to have a potential relevance for the physics of protoneu- 
tron stars where neutrinos are trapped during the first few seconds of the 
stellar evolution. 

The effect of neutrino trapping on color-superconducting quark matter 
has been previously discussed in Ref. |3| There it was found that a nonzero 
neutrino chemical potential favors the 2SC phase and disfavors the CFL 
phase. This is not unexpected because the neutrino chemical potential is 
related to the conserved lepton number in the system and therefore it also 
favors the presence of (negatively) charged leptons. This helps 2SC-type 
pairing because electrical neutrality in quark matter can be achieved with- 
out inducing a very large mismatch between the Fermi surfaces of up and 
down quarks. The CFL phase, on the other hand, i s el ectrically and color 
neutral in the absence of charged leptons when T = A nonzero neutrino 
chemical potential can only spoil CFL-type pairing. 

A more systematic survey of the phase diagram in the space of tem- 
perature, quark and lepton-number chemical potentials was performed in 
Ref. El In particular, this included the possibility of gapless phases, which 
have not been taken into account in Ref. 01 While such phases are gener- 
ally unstable at zero temperature this is not always the case at nonzero 
temperature!^^ Keeping this in mind, we shall merely localize the "prob- 
lematic" regions in the phase diagram, where unconventional pairing is 
unavoidable. We shall refrain, however, from speculating on various pos- 
sibilities for the true ground state (see, e.g., Refs. OUlini 1^ 1^ 1771 1^ . 
since these are still under debate. 

In the application to protoneutron stars, it is of interest to cover a 
range of parameters that could provide a total lepton fraction in quark 
matter of up to about 0.4. This is the value of the lepton-to-baryon charge 
ratio in iron cores of progenitor stars. Because of the conservation of both 
lepton and baryon charges, this value is also close to the lepton fraction in 
protoneutron stars at early times, when the leptons did not have a chance 
to diffuse through dense matter and escape from the star.^ In the model 
introduced in the next section, almost the whole range of possibilities will 
be covered by restricting the quark chemical potential to /i < 500 MeV and 
the neutrino chemical potential to fi^^ < 400 MeV. 
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4 Riister, Werth, Buballa, Shovkovy, and Rischke 

2. Model 

Let us start by introducing the effective model of QCD used in the analysis. 
This is a three-flavor quark model with a local NJL-type interaction, whose 
Lagrangian density is given by 

8 

- K jdet [i> (1 + 75) ^] + det (1 - 75) | , (1) 

where the quark spinor field ip"^ carries color (a — r,g,b) and flavor 
(a — u, d, s) indices. The matrix of quark current masses is given by 
TO — diag^(TO„, TOd, TOs). Regarding other notations, Aa with a = 1, ... ,8 
are the Gell-Mann matrices in flavor space, and Aq = \/2/3 1/. The charge 
conjugate spinors are defined as follows: = CV^^ and V'c = V'^C'? where 
^ = i^t^o jg -j-jj^g Dirac conjugate spinor and C = 17^7" is the charge conju- 
gation matrix. 

The model in Eq. ^ should be viewed as an effective model of strongly 
interacting matter that captures at least some key features of QCD dynam- 
ics. The Lagrangian density contains three different interaction terms which 
are chosen to respect the symmetries of QCD. Note that we include the 't 
Hooft interaction whose strength is determined by the coupling constant 
K. This term breaks the U{\) axial symmetry. 

The term in the second line of Eq. describes a scalar diquark inter- 
action in the color antitriplet and flavor antitriplet channel. For symmetry 
reasons there should also be a pseudoscalar diquark interaction with the 
same coupling constant. This term would b e im portant to describe Gold- 
stone boson condensation in the CFL phase.EHI 

We use the following set of model parameters!^ 

m„ = TOd 5.5 MeV, m, = 140.7 MeV, (2a) 
GsA^ = 1.835, KK^ = 12.36, A = 602.3 MeV. (2b) 

After fixing the masses of the up and down quarks at equal values, mu.d = 
5.5 MeV, the other four parameters are chosen to reproduce the following 
four observables of vacuum QCD: to^ = 135.0 MeV, rriK = 497.7 MeV, 
TO,,' = 957.8 MeV, and = 92.4 MeV.ElThis parameter set gives = 
514.8 MeV. 
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Phase diagram of neutral quark matter at moderate densities 5 

In Ref.|^ the diquark couphng Go was not fixed by the fit of the meson 
spectrum in vacuum. In general, it is expected to be of the same order as 
the quark- antiquark coupling Gs- Here, we choose the coupling strength 
with Gd = jGs which follows from the Fierz identity. 

The grand partition function, up to an irrelevant normalization con- 
stant, is given by 

Z = e-^V/T^ J p^p^ ^^ J, (£+^/i7 V) ^ (3) 

where fl is the thermodynamic potential density, V is the volume of the 
three-space, and /t is a diagonal matrix of quark chemical potentials. In 
chemical equilibrium (which provides f3 equilibrium as a special case), the 
nontrivial components of this matrix are extracted from the following rela- 
tion: 

= (m^-^ + MQQf ) Sat + [m3 (Ta),, + iTs)J S'^" ■ (4) 

Here fi is the quark chemical potential (by definition, /i = ij-b/3 where fis 
is the baryon chemical potential), /ig is the chemical potential of electric 
charge, while and ns are color chemical potentials associated with two 
mutually commuting color charges of the SU{3)c gauge group, cf. Ref. |32 
The explicit form of the electric charge matrix is Qf = diag^(|, — i, — i), 
and the explicit form of the color charge matrices is T3 = diag^(i, — i, 0) 
and\/3r8 = diag,(i,i,-l). 

In order to calculate the mean-field thermodynamic potential at tem- 
perature T, we first linearize the interaction in the presence of the diquark 
condensates Ac ~ {i^c)a^l5^°'^'^^abci^'i3 (no sum over c) and the quark- 
antiquark condensates a a ^ V-'aV'S (no sum over a). Then, integrating out 
the quark fields and neglecting the fiuctuations of composite order parame- 
ters, we arrive at the following expression for the thermodynamic potential: 

3 3 

c— 1 a— 1 

IV ^ 1 

where we also added the contribution of Icptons, fi^, which will be specified 
later. 

We should note that we have restricted ourselves to field contractions 
corresponding to the Hartree approximation. In a more complete treatment. 
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6 Riister, Werth, Buballa, Shovkovy, and Rischke 

among others, the 't Hooft interaction term gives also rise to mixed con- 
tributions contai ning both diquark and quark- antiquark condensates, i.e., 
cx ''"al^aP-^^ In this study, we neglect such terms for simplicity. 

While their presence may change the results quantitatively, one does not 
expect them to modify the qualitative structure of the phase diagram. 

In Eq. lO, S*^^ is the inverse full quark propagator in the Nambu-Gorkov 
representation. 



with the diagonal elements being the inverse Dirac propagators of quarks 
and of charge- conjugate quarks, 

[G±]-i-7^^M±/'7o-M, (7) 

where = (fco,k) denotes the four-momentum of the quark. At nonzero 
temperature, we use the Matsubara imaginary time formalism. Therefore, 
the energy fco is replaced with —iujn where u;„ = (2n+l)7rT are the fermionic 
Matsubara frequencies. Accordingly, the sum over K in Eq. |SJ) should 
be interpreted as a sum over integer n and an integral over the three- 
momentum k. 

The constituent quark mass matrix M — diagj(M„, M^, Ms) is defined 

by 

Ma — rria — 4:Gs<Ja + 2Kapa^ , (8) 

where are the quark- antiquark condensates, and the set of indices 
(a, /?, 7) is a permutation of (u, d, s). 

The off-diagonal components of the inverse propagator © are the so- 
called gap matrices given in terms of three diquark condensates. The color- 
flavor structure of these matrices is given by 

(*-):f=-E^"''^^-^=75, (9) 

c 

and = 7°($~)^7°. Here, as before, a and b refer to the color components 
and a and P refer to the flavor components. Hence, the gap parameters Ai, 
A2, and A3 correspond to the down-strange, the up-strange and the up- 
down diquark condensates, respectively. All three of them originate from the 
color-antitriplet, flavor-antitriplet diquark pairing channel. For simplicity, 
the color and flavor symmetric condensates are neglected in this study. They 
were shown to be small and not crucial for the qualitative understanding 
of the phase diagram.^^ 
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The determinant of the inverse quark propagator can be decomposed as 
followsPl 

1=1 ^ ^ 

where are eighteen independent positive energy eigenvalues. The Mat- 
subara su mma tion in Eq. |SJ| can then be done analytically by employing 
the relatioi]^ 

n ^ ^ 

Then, we arrive at the following mean-field expression for the pressure: 
P=-^=2^Y. rffcfc' [|e.| +2rin(l + e-^) 

, 3 3 
c—1 a—1 

T ^ ^ r / B,-^. 



5EE / dkkHn{l 



+ + . (12) 

where the contributions of electrons and muons with masses rrie « 0.511 
MeV and « 105.66 MeV, as well as the contributions of neutrinos 
were included. Note that muons may exist in matter in /3 equilibrium and, 
therefore, they are included in the model for consistency. However, being 
about 200 times heavier than electrons, they do not play a big role in the 
analysis. 

3. Phase diagram in absence of neutrino trapping 

In this section, we consider the case without neutrino trapping in quark 
matter. This is expected to be a good approximation for matter inside a 
neutron star after the short deleptonization period is over. 

The expression for the pressure in Eq. H12() has a physical meaning only 
when the chiral and color-superconducting order parameters, CTq, and Ac, 
satisfy the following set of six gap equations: 
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8 Riister, Werth, Buballa, Shovkovy, and Rischke 

To enforce the conditions of local charge neutrality in dense matter, we also 
require three other equations to be satisfied, 

dp ^ dp ^ dp ^ 

uq = - — = , na = - — = , 718 = ^ — = . (14) 
d^iQ d^i^ dfis 

These fix the values of the three corresponding chemical potentials, /ig, ^3 
and /^8- After these are fixed, only the quark chemical potential /i is left as 
a free parameter. 

In order to obtain the phase diagram, one has to find the ground state 
of matter for each given set of the parameters in the model. In the case of 
locally neutral matter, there are two parameters that should be specified: 
temperature T and quark chemical potential /z. After these are fixed, one 
has to compare the values of the pressure in all competing neutral phases of 
quark matter. The ground state corresponds to the phase with the highest 
pressure. 

By using standard numerical recipes, it is not extremely difficult to find 
a solution to the given set of nine nonlinear equations. Complications arise, 
however, due to the fact that often the solution is not unique. The existence 
of different solutions to the same set of equations, (|13|l and (|14|l . reflects the 
physical fact that there could exist several competing neutral phases with 
different physical properties. Among these phases, all but one are unstable 
or metastable. In order to take this into account, one should look for the 
solutions of the following 8 types: 

(1) Normal quark (NQ) phase: Ai = A2 = A3 = 0; 

(2) 2SC phase: Ai = A2 = 0, and A3 ^ 0; 

(3) 2SCus phase: Ai = A3 = 0, and A2 ^ 0; 

(4) 2SCds phase: A2 = A3 = 0, and Ai ^ 0; 

(5) uSC phase: A2 ^ 0, A3 ^ 0, and Ai = 0; 

(6) dSC phase: Ai ^ 0, A3 ^ 0, and A2 = 0; 

(7) sSC phase: Ai 7^ 0, A2 7^ 0, and A3 = 0; 

(8) CFL phase: Ai ^ 0, A2 ^ 0, A3 ^ 0. 

Then, we calculate the values of the pressure in all nonequivalent phases, 
and determine the ground state as the phase with the highest pressure. After 
this is done, we study additionally the spectrum of low-energy quasiparticles 
in search for the existence of gapless modes. This allows us to refine the 
specific nature of the ground state. 

In the above definition of the eight phases in terms of Ac, we have 
ignored the quark-antiquark condensates CTq.. In fact, in the chiral limit 
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Phase diagram of neutral quark matter at moderate densities 9 

(nia — 0), the quantities cTq, are good order parameters and we could define 
additional sub-phases characterized by nonvanishing values of one or more 
CTq,. With the model parameters at hand, however, chiral symmetry is broken 
explicitly by the nonzero current quark masses, and the values of CTq never 
vanish. Hence, in a strict sense it is impossible to define any new phases in 
terms of a a ■ 

Of course, this does not exclude the possibility of discontinuous changes 
in (Tq, at some line in the plane of temperature and quark chemical poten- 
tial, thereby constituting a first-order phase transition line. It is generally 
expected that the "would-be" chiral phase transition remains first-order 
at low temperatures, even for nonzero quark masses. Above some critical 
temperature, however, this line could end in a critical endpoint and there is 
only a smooth crossover at higher tem peratures. This picture emerges from 
NJL-model studies, both, withoulPSI and withPSI diquark pairing. 

The numerical results for neutral quark matter are summarized in Fig.^ 
This shows the phase diagram in the plane of temperature and quark chem- 
ical potential, obtained in the mean-field approximation in model (0 . The 
corresponding dynamical quark masses, gap parameters, and three charge 
chemical potentials are displayed in Fig. [21 All quantities are plotted as func- 
tions of fi for three different values of the temperature: T = 0, 20, 40 MeV. 

In the region of small quark chemical potentials and low temperatures, 
the phase diagram is dominated by the normal phase in which the ap- 
proximate chiral symmetry is broken, and in which quarks have relatively 
large constituent masses. This is denoted by xSB in Fig.Q With increasing 
the temperature, this phase changes smoothly into the NQ phase in which 
quark masses are relatively small. Because of explicit breaking of the chiral 
symmetry in the model at hand, there is no need for a phase transition 
between the two regimes. 

However, as pointed out above, the symmetry argument does not ex- 
clude the possibility of an "accidental" (first-order) chiral phase transition. 
As expected, at lower temperatures we find a line of first-order chiral phase 
transitions. It is located within a relatively narrow window of the quark 
chemical potentials (336 MeV ^ M ^ 368 MeV) which are of the order 
of the vacuum values of the light-quark constituent masses. (For the pa- 
rameters used in the calculation one obtains Af„ — Md — 367.7 MeV and 
Ms = 549.5 MeV in vacuum!--!) At this critical line, the quark chiral con- 
densates, as well as the quark constituent masses, change discontinuously. 
With increasing temperature, the size of the discontinuity decreases, and 
the line terminates at the endpoint located at {Tcr,fJ,ci) ~ (56,336) MeV, 
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Fig. 1. The phase diagram of neutral quark matter in model First-order phase 
boundaries are indicated by bold solid lines, whereas the thin solid lines mark second- 
order phase boundaries between two phases which differ by one or more nonzero diquark 
condensates. The dashed lines indicate the (dis-)appearance of gapless modes in different 
phases, and they do not correspond to phase transitions. 



see Fig. ^ 

The location of the critical endpoint is consistent with other mean-field 
studies of NJL models with similar sets of parameters. "^"^ "^^ This agreement 
does not need to be exact because, in contrast to the studies in Refs. 1351 
1361 here we imposed the condition of electric charge neutrality in quark 
matter. (Note that the color neutrality is satisfied automatically in the 
normal phase.) One may argue, however, that the additional constraint of 
neutrality is unlikely to play a big role in the vicinity of the endpoint. 

It is appropriate to mention here that the location of the critical end- 
point might be affected very much by fluctuations of the composite chiral 
fields. These are not included in the mean-field studies of the NJL model. 
In fact, this is probably the main reason for their inability to pin down 
the locatio n of the critical endpoint consistent, for example, with lattice 
calculations!^ (It is fair to mention that the current lattice calculations 
are not very reliable at nonzero ^ either.) Therefore, the predictions of this 
study, as well as of those in Refs. 1351 and Kffil regarding the critical endpoint 
cannot be considered as very reliable. 

When the quark chemical potential exceeds some critical value and the 
temperature is not too large, a Cooper instability with respect to diquark 
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Fig. 2. Dependence of the quark masses, of the gap parameters, and of the electric and 
color charge chemical potentials on the quark chemical potential at a fixed temperature, 
T = MeV (three upper panels), T = 20 MeV (three middle panels), and T = 40 MeV 
(three lower panels). 



condensation should develop in the system. Without enforcing neutrality, 
i.e., if the chemical potentials of up and down quarks are equal, this happens 
immediately after the chiral phase transition when the density becomes 
nonzero.!^ In the present model, this is not the case at low temperatures. 
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In order to understand this, one should inspect the various quantities 
at T = which are displayed in the upper three panels of Fig. [3 At the 
chiral phase boundary, the up and down quark masses become relatively 
small, whereas the strange quark mass experiences only a moderate drop 
of about 84 MeV induced by the 't Hooft interaction. This is not sufficient 
to populate any strange quark states at the given chemical potential, and 
the system mainly consists of up and down quarks together with a small 
fraction of electrons, see Fig.|3| The electric charge chemical potential which 
is needed to maintain neutrality in this regime is between about —73 and 
—94 MeV. It turns out that the resulting splitting of the up and down 
quark Fermi momenta is too large for the given diquark coupling strength 
to enable diquark pairing and the system stays in the normal phase. 

10 



1 



0.01 



0.001 

300 350 400 450 500 

\i [MeV] 

Fig. 3. The dependence of the number densities of quarks and electrons on the quark 
chemical potential at T = MeV. Note that the densities of all three quark flavors 
coincide above fi = 457 MeV. The density of muons vanishes for all values of /i. 

At /i w 432 MeV, the chemical potential felt by the strange quarks, 
M~Mq/3, reaches the strange quark mass and the density of strange quarks 
becomes nonzero. At first, this density is too small to play a sizable role in 
neutralizing matter, or in enabling strangc-nonstrange cross-flavor diquark 
pairing, see Fig. The NQ phase becomes metastable against the gapless 
CFL (gCFL) phase at /LtgCFL ~ 443 MeV. This is the point of a first-order 
phase transition. It is marked by a drop of the strange quark mass by 
about 121 MeV. As a consequence, strange quarks become more abundant 
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Phase diagram of neutral quark matter at moderate densities 13 

and pairing gets easier. Yet, in the gCFL phase, the strange quark mass 
is still relatively large, and the standard BCS pairing between strange and 
light (i.e., up and down) quarks is not possible. In contrast to the regular 
CFL phase, the gCFL phase requires a nonzero density of electrons to stay 
electrically neutral. At T = 0, therefore, one could use the value of the 
electron density as a formal order parameter that distinguishes these two 
phases .1^ 

With increasing the chemical potential further (still at T = 0), the 
strange quark mass decreases and the cross-flavor Cooper pairing gets 
stronger. Thus, the gCFL phase eventually turns into the regular CFL 
phase at hcfl ~ 457 MeV. The electron density goes to zero at this point, 
as it should. This is indicated by the vanishing value of fig in the CFL 
phase, see the upper right panel in Fig. [3 We remind that the CFL phase 
is neutral because of having equal number densities of all three quark fla- 
vors, Tiu ^ rid = Us, see Fig. O This equality is enforced by the pairing 
mech anis m, and this is true even when the quark masses are not exactly 
equal P21 

Let us mention here that the same NJL model at zero temperature was 
studied previously in Ref.'S' The results of Ref.|2] agree with those presented 
here only when the quark chemical potential is larger than the critical value 
for the transition to the CFL phase at 457 MeV. The appearance of the 
gCFL phase for 443 A* 457 MeV was not recognized in Ref.'B', however. 
Instead, it was suggested that there exists a narrow (about 12 MeV wide) 
window of values of the quark chemical potential around fj, « 450 MeV in 
which the 2SC phase is the ground state. By examining the same region, 
we find that the 2SC phase does not appear there. 

This is illustrated in Fig.0]where the pressure of three different solutions 
is displayed. Had we ignored the gCFL solution (thin solid line), the 2SC 
solution (dashed line) would indeed be the most favored one in the interval 
between fi « 445 MeV and « 457 MeV. After including the gCFL phase 
in the analysis, this is no longer the case. 

Now let us turn to the case of nonzero temperature. One might suggest 
that this should be analogous to the zero-temperature case, except that 
Cooper pairing is somewhat suppressed by thermal effects. In contrast to 
this naive expectation, the thermal distributions of quasiparticles together 
with the local neutrality conditions open qualitatively new possibilities that 
were absent at T = 0. As in the case of the two-flavor model of Ref. |S1 a 
moderate thermal smearing of mismatched Fermi surfaces could increase 
the probability of creating zero-momentum Cooper pairs without running 
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Fig. 4. Pressure divided by /i'* for different neutral solutions of the gap equations at 
T = as functions of the quark chemical potential fi: regular CFL (bold solid line), 
gapless CFL (thin solid line), 2SC (dashed line), normal (dotted line). 

into a conflict with Pauli blocking. This leads to the appearance of several 
color-superconducting phases that could not exist at zero temperatures. 

With increasing the temperature, the first qualitatively new feature in 
the phase diagram appears when 5 < T < 10 MeV. In this temperature in- 
terval, the NQ phase is replaced by the uSC phase when the quark chemical 
potential exceeds the critical value of about 444 MeV. The corresponding 
transition is a first-order phase transition, see Fig.^ Increasing the chemi- 
cal potential further by several MeV, the uSC phase is then replaced by the 
gCFL phase, and the gCFL phase later turns gradually into the (m)CFL 
phase. (In this study, we do not distinguish between the CFL phase and the 
mCFL phase.l^) Note that, in the model at hand, the transition between 
the uSC and the gCFL phase is of second order in the following two tem- 
perature intervals: 5 < T < 9 MeV and T >24 MeV. On the other hand, it 
is a first-order transition when 9 < T < 24 MeV. Leaving aside its unusual 
appearance, this is likely to be an "accidental" property in the model for a 
given set of parameters. 

The tra nsitio n from the gCFL to the CFL phase is a smooth crossover 
at all T oW^^The reason is that the electron density is not a good order 
parameter that could be used to distinguish the gCFL from the CFL phase 
when the temperature is nonzero. This is also confirmed by the numerical 
results for the electric charge chemical potential /xq in Fig. El While at zero 
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temperature the value of /ig vanishes identically in the CFL phase, this is 
not the case at nonzero temperatures. 

Another new feature in the phase diagram appears when the tempera- 
ture is above about 11 MeV. In this case, with increasing the quark chem- 
ical potential, the Cooper instability happens immediately after the x^B 
phase. The corresponding critical value of the quark chemical potential is 
rather low, about 365 MeV. The first color-superconducting phase is the 
gapless 2SC (g2SC) phased This phase is replaced with the 2SC phase in 
a crossover transition only when /i > 445 MeV. The 2SC is then followed 
by the gapless uSC (guSC) phase, by the uSC phase, by the gCFL phase 
and, eventually, by the CFL phase (see Fig.P). 

In the NJL model at hand, determined by the parameters in Eq. ||2Jl, 
we do not find the dSC phase as the ground state anywhere in the phase 
diagram. This is similar to the conclusion of Refs. [3 and I14L but differs 
from that of Refs. '101 andll3l This should not be surprising because, as was 
noted earlier, -'^^ -'^^Jthe appearance of the dSC phase is rather sensitive to 
a specific choice of parameters in the NJL model. 

The phase diagram in Fig.QJhas a very specific ordering of quark phases. 
One might ask if this ordering is robust against the modification of the 
parameters of the model at hand. We can argue that some features are 
indeed quite robust, while others are not."'-^ 

It should be clear that the appearance of color-superconducting phases 
under the stress of neutrality constraints is very sensitive to the strength of 
diquark coupling. In the case of two-flavor quark matter, this was demon- 
strated very clearly in Rcf. El at zero as well as at nonzero te mpera tures. 
The same statement remains true in three-flavor quark matter! ' ■ ' I '^ l 

4. Phase diagram in presence of neutrino trapping 

In the case of neutrino trapping, the chemical potentials of individual quark 
and lepton species can be expressed in terms of six chemical potentials ac- 
cording to their content of conserved charges. For the quarks, which carry 
quark number, color and electric charge, the chemical potentials were in- 
troduced in Eq. For the neutrinos, which carry only lepton number, 
the chemical potentials are 

M.^, = /^Le : Mi^^, = fJ-L^ ■ (15) 

while for the electrons and muons, carrying both lepton number and electric 
charge, the chemical potentials read 
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In order to obtain the phase diagram, one has to determine the ground state 
of matter for each given set of the parameters. In the case of locally neutral 
matter with trapped neutrinos, there are four parameters that should be 
specified: the temperature T, the quark chemical potential /i as well as the 
two lepton family chemical potentials (II^ and fiL^ ■ After these are fixed, 
the values of the pressure in all competing neutral phases of quark matter 
should be compared. The phase with the largest pressure is the ground 
state. For our purposes, it is sufficient to take the vanishing muon lepton- 
number chemical potential, i.e., fj,L^ — 0. This is expected to be a good 
approximation for matter inside a protoncutron star. 

4.1. General effect of neutrino trapping 

As mentioned in the Introduction, neutrino trapping favors the 2SC phase 
and disfavors the CFL phase!^ This is a consequence of the modified /3- 
equilibrium condition in the system. In this section, we would like to em- 
phasize that this is a model- independent effect. In order to understand the 
physics behind it, it is instructive to start our consideration from a very 
simple toy model. Many of its qualitative features are also observed in a 
self-consistent numerical analysis of the NJL model. 

Let us first assume that strange quarks are very heavy and consider a 
gas of non-interacting massless up and down quarks in the normal phase at 
T = 0. As required by /3 equilibrium, electrons and electron neutrinos are 
also present in the system. (Note that in this section we neglect muons and 
muon neutrinos for simplicity.) 

In the absence of Cooper pairing, the densities of quarks and leptons 
are given by 

nu.d = , ne = — ^ , Ui,^ = — . (17) 

Expressing the chemical potentials through /i, /ig and (J-l^, and imposing 
electric charge neutrality, one arrives at the following relation: 

2(1 + ^y)' - (i-^y)' -{x-yf ^ 0, (18) 

where we have introduced the chemical potential ratios x — and 
y — fiQ / II. The above cubic equation can be solved for y (electric chemical 
potential) at any given x (lepton- number chemical potential). The result 
can be used to calculate the ratio of quark chemical potentials, iJ-d/fJ-u = 
(3-y)/(3 + 22/). 
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The ratio fid/fJ-u as a function of fJ-L^/l^ is shown in Fig.|Sl At vanishing 
/i/,^, one finds y « —0.219 and, thus, fid/fJ-u ~ 1.256 (note that this value is 
very close to 2^/^ « 1.260). This result corresponds to the following ratios 
of the number densities in the system: Uu/rid w 0.504 and Ue/nd ~ 0.003, 
reflecting that the density of electrons is tiny and the charge of the up 
quarks has to be balanced by approximately twice as many down quarks. 

1.3 
1.25 

1.2 
1.15 

=1 

1.1 

1.05 
1 

0.2 0.4 0.6 0.8 1 

Fig. 5. Ratio of down and up quark chemical potentials as a function of ^l^/a* in the 
toy model. The crosses mark the solutions at several values of the lepton fraction. 

At /Xi^ — II, on the other hand, the real solution to Eq. H18() is y = 0, 
i.e., the up and down Fermi momenta become equal. This can be seen most 
easily if one inverts the problem and solves Eq. (|18|l for x at given y. When 
y = one finds x = I, meaning that fid — fJ'u and, in turn, suggesting that 
pairing between up and down quarks is unobstructed at fiL^ — This is 
in contrast to the case of vanishing fiL^, when the two Fermi surfaces are 
split by about 25%, and pairing is difficult. 

It is appropriate to mention that many features of the above consid- 
erations would not change much even when Cooper pairing is taken into 
account. The reason is that the corresponding corrections to the quark 
densities are parametrically suppressed by a factor of order (A/^)^. 

In order to estimate the magnitude of the effect in the case of dense 
matter in protoneutron stars, we indicate several typical values of the lepton 
fractions Y^^ in Fig. [3 As mentioned earlier, Yl^ is expected to be of order 
0.4 right after the collapse of the iron core of a progenitor star. According to 
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Fig- El this corresponds to /id/Mu ~ 1-1: while the spHtting between the 
up and down Fermi surfaces does not disappear completely, it gets reduced 
considerably compared to its value in the absence of trapped neutrinos. 
This reduction substantially facilitates the cross-flavor pairing of up and 
down quarks. The effect is gradually washed out during about a dozen of 
seconds of the deleptonization period when the value of Yl^ decreases to 
zero. 

The toy model is easily modified to the opposite extreme of three mass- 
less quark flavors. Basically, this corresponds to replacing Eq. H18|) by 



In the absence of neutrino trapping, x = 0, the only real solution to this 
equation is ?; 0, indicating that the chemical potentials (which also coin- 
cide with the Fermi momenta) of up, down, and strange quarks are equal. 
This reflects the fact that the system with equal densities of up, down, 
and strange quarks is neutral by itself, without electrons. With increasing 
X (X fJ.L^, the solution requires a nonzero y cx /iq, suggesting that up-down 
and up-strange pairing becomes more difficult. To see this more clearly, we 
can go one step further in the analysis of the toy model. 

Let us assume that the quarks are paired in a regular, i.e., fully gapped, 
CFL phase at T = 0. Then, as shown in Ref. the quark part of the 
matter is automatically electrically neutral. Hence, if we want to keep the 
whole system electrically and color neutral, there must be no electrons. 
Obviously, this is easily realized without trapped neutrinos by setting /ig 
equal to zero. At non- vanishing fj.^^ the situation is more complicated. The 
quark part is still neutral by itself and therefore no electrons are admitted. 
Hence, the electron chemical potential fj,e = ^il^ — must vanish, and 
consequently /xq should be nonzero and equal to /iL^. It is natural to ask 
what should be the values of the color chemical potentials /13 and /ig in the 
CFL phase when 7^ 0. 

In order to analyze the stress on the CFL phase due to nonzero /i^^, 
we follow the same approach as in Refs. El and El In this analytical con- 
sideration, we also account for the effect of the strange quark mass simply 
by shifting the strange quark chemical potential by — M^/(2/z). In our no- 
tation, pairing of CFL-type requires the following "common" values of the 



2{l +'^yf 2{l -^-yf {x - yf 
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(19) 
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Fermi momenta of paired quarks: 



2 



M' 

common _ s /'OHq^ 

PF,(ru,gd,bs) " M - : {Ma.) 

common _ ,, i MQ i ^8 (0(\W\ 
PF,{rd,gu) ~ ^ ^ 

common _ , MQ , M3 MS 

common _ ,, MQ M3 M8 ^r,r,j\ 

These are used to calculate the pressure in the toy model, 

7r2 127r2 247r2 ^ ' 

By making use of this expression, one easily derives the neutrality conditions 
as in Eq. (|14|l . In order to solve them, it is useful to note that 

nQ-n,-—ns = . (22) 

Thus, it becomes obvious that charge neutrality requires /xg = ^ll^- The 
other useful observation is that the expression for is proportional to 
M3 + MQ- So, it is vanishing if (and only if) = — /^q, which means that 
/13 = —fJ.L^ ■ Finally, one can check that the third neutrality condition ng — 
requires 

The results for the charge chemical potentials /^g, and /xg imply the 
following magnitude of stress on pairing in the CFL phase: 

SlJ.(rd,gu) = ^ ^'^ = Miei (24a) 

Sf^irsM) = ^^^y^ = + (24b) 

'5M(s.,6d) = = (24c) 

Note that there is no mismatch between the values of the chemical poten- 
tials of the other three quarks, fi"!^ — fig — fj.^ — ^ ~ M^/{6fi). 
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From Eq. (|24|l we see that the largest mismatch occurs in the (rs, bu) 
pair (for positive ^l^)- The CFL phase can withstand the stress only if 
the value of (5/i(rs,bu) is less than A2. A larger mismatch should drive a 
transition to a gapless phase exactly as in Refs. Eland El Thus, the critical 
value of the lepton-number chemical potential is 



When iiL^ > IJ, " , the CFL phase turns into the gCFL' phase, which is a 
variant of the gCFL phased By definition, the gapless mode with a linear 
dispersion relation in the gCFL' phase is rs-bu instead of gs-bd as in the 
standard gCFL phase. (Let us remind that the mode aa-b(3 is defined by its 
dispersion relation which interpolates between the dispersion relations of 
hole-type excitations of aa-quark at small momenta, fc ^ , and particle- 
type excitations of 6/3-quark at large momenta, fc 3> /if.) 

In order to see what this means for the physics of protoneutron stars, we 
should again try to relate the value of to the lepton fraction. As we have 
seen, there are no electrons in the (regular) CFL phase at T = 0. Therefore, 
the entire lepton number is carried by neutrinos. For the baryon density 
we may neglect the pairing effects to first approximation and employ the 
ideal-gas relations. This yields 



Inserting typical numbers, /j, > 400 MeV and ^j^^ ^ A < 100 MeV, one finds 
Yl^ < 10^"^. Thus, there is practically no chance to find a sizable amount of 
leptons in the CFL phase. The constraint gets relaxed slightly at nonzero 
temperatures and/or in the gCFL phase, but the lepton fraction remains 
rather small even then (our numerical results indicate that, in general, 
Yl, < 0.05 in the CFL phase). 

4.2. Three-dimensional phase diagram 

The simple toy-model considerations in the previous subsection give a qual- 
itative understanding of the effect of neutrino trapping on the mismatch 
of the quark Fermi momenta and, thus, on the pairing properties of two- 
and three-flavor quark matter. Now, we turn to a more detailed numerical 
analysis of the phase diagram in the NJL model defined in Sec. El 

The general features of the phase diagram in the three-dimensional 
space, spanned by the quark chemical potential /U, the lepton-number chem- 




(26) 
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ical potential /iL^: the temperature T, are depicted in Fig.|^ Because 
of a rather comphcated structure of the diagram, only the four main phases 
(xSB, NQ, 2SC and CFL) are shown explicitly. Although it is not labeled, 
a thin slice of a fifth phase, the uSC phase, squeezed in between the 2SC 
and CFL phases, can also be seen. While lacking detailed information, the 
phase diagram in Fig. El gives a clear overall picture. Among other things, 
one sees, for example, that the CFL phase becomes strongly disfavored with 
increasing fi^^ and gets gradually replaced by the 2SC phase. 



/iL, [MeV] - - - 



100 



T [MeV 




400 
IJ [MeV] 



450 



500 



Fig. 6. General structure of the phase diagram of neutral dense quark matter in the 
three-dimensional space spanned by the quark chemical potential /i, the lepton-number 
chemical potential fij^ , and the temperature T. 



In order to discuss the structure of the phase diagram in more detail we 
proceed by showing several two-dimensional slices of it. These are obtained 
by keeping one of the chemical potentials, fi or /il^, fixed and varying the 
other two parameters. 
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4.3. T fi phase diagram 

We begin with presenting the phase diagrams at two fixed values of the 
lepton- number chemical potential, /ii^ = 200 MeV (upper panel) and fJ.L^ = 
400 MeV (lower panel), in Fig.[71 The general effects of neutrino trapping 
can be understood by analyzing the similarities and differences between the 
diagrams with and without neutrino trapping. 
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Fig. 7. The phase diagrams of neutral quark matter at fixed lepton-number chemical 
potentials fj.^^ = 200 MeV (upper panel), and /i^^ = 400 MeV (lower panel), cf. Fig. Ill 
in the absence of neutrino trapping. 
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Here it is appropriate to note that a schematic version of the T-fi phase 
diagram at fiL^ = 200 MeV was first presented in Ref. El see the right 
panel of Fig. 4 there. If one ignores the complications due to the presence 
of the uSC phase and various gapless phases, the results of Ref. El are in 
qualitative agreement with the results presented here. 

In order to understand the basic characteristics of different phases in the 
phase diagrams in Fig.[7| we also present the results for the dynamical quark 
masses, the gap parameters, and the charge chemical potentials. These are 
plotted as functions of the quark chemical potential in Figs. |H1 and El for 
two different values of the temperature in the case of — 200 MeV and 
Hl^ — 400 MeV, respectively. 
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Fig. 8. Dependence of the quark masses, of the gap parameters, and of the electric and 
color charge chemical potentials on the quark chemical potential at a fixed temperature, 
T = MeV (three upper panels) and T = 40 MeV (three lower panels). The lepton- 
number chemical potential is kept fixed a,t fij^^ = 200 MeV. 



As in the absence of neutrino trapping, see Fig. ^ there are roughly four 
distinct regimes in the phase diagrams in Fig. |3 At low temperature and 
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Fig. 9. Dependence of the quark masses, of the gap parameters, and of the electric and 
color charge chemical potentials on the quark chemical potential at a fixed temperature, 
T = MeV (three upper panels), and T = 40 MeV (three lower panels). The lepton- 
number chemical potential is kept fixed at fj,j^^ = 400 MeV. 



low quark chemical potential, there is the xSB region. Here quarks have 
relatively large constituent masses which are close to the vacuum values, 
see Figs, island El We note that the xSB phase is rather insensitive to the 
presence of a nonzero lepton-number chemical potential. With increasing 
Hi^^ the phase boundary is only slightly shifted to lower values of fj,. This 
is just another manifestation of the strengthening of the 2SC phase due to 
neutrino trapping. 

With increasing temperature, the x^B phase turns into the normal 
quark matter phase, whose qualitative features are little affected by the 
lepton-number chemical potential. 

The third regime is located at relatively low temperatures but at quark 
chemical potentials higher than in the x^B phase. In this region, the masses 
of the up and down quarks have already dropped to values well below their 
respective chemical potentials while the strange quark mass is still large, 
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see left columns of panels in Figs. |H1 and |^ As a consequence, up and down 
quarks are quite abundant but strange quarks are essentially absent. 

It turns out that the detailed phase structure in this region is very sen- 
sitive to the lepton-number chemical potential. At — 0, the pairing 
between up and down quarks is strongly hampered by the constraints of 
neutrality and f3 equilibrium, see Fig. ^ The situation changes dramati- 
cally with increasing the value of the lepton-number chemical potential, see 
Fig. \7\ The low-temperature region of the normal phase of quark matter 
is replaced by the (g)2SC phase (e.g., at ^ = 400 MeV, this happens at 
fiL^ ~ 110 MeV). With /i^^ increasing further, no qualitative changes hap- 
pen in this part of the phase diagram, except that the area of the (g)2SC 
phase expands slightly. 

Finally, the region in the phase diagram at low temperatures and large 
quark chemical potentials corresponds to phases in which the cross-flavor 
strange-nonstrange Cooper pairing becomes possible. In general, as the 
strength of pairing increases with the quark chemical potential, the sys- 
tem passes through regions of the gapless uSC (guSC), uSC, and gCFL 
phases and finally reaches the CFL phase. (Of course, the intermediate 
phases may not always be realized.) The effect of neutrino trapping, which 
grows with increasing the lepton-number chemical potential, is to push out 
the location of the strange-nonstrange pairing region to larger values of /i. 
Of course, this is in agreement with the general arguments in Sec. 14.11 

As we mentioned earlier, the presence of the lepton-number chemical 
potential /i^^ leads to a change of the quark Fermi momenta. This change 
in turn affects Cooper pairing of quarks, facilitating the appearance of some 
phases and suppressing others. As it turns out, there is also another qual- 
itative effect due to a nonzero value of fj,L^- In particular, we find several 
new variants of gapless phases which did not exist at vanishing /iL^- In 
Figs. 171 and im these are denoted by the same names, g2SC or gCFL, but 
with one or two primes added. 

We define the g2SC' as the gapless two-flavor color-superconducting 
phase in which the gapless excitations correspond to rd-gu and gd-ru 
modes instead of the usual ru~gd and gu~rd ones, i.e., u and d flavors are 
exchanged as compared to the usual g2SC phase. The g2SC' phase becomes 
possible only when the value of (/i™ — /igd)/2 = (/ig + ^,^)/2 is positive 
and larger than A3. The other phases are defined in a similar manner. In 
particular, the gCFL' phase, already introduced in Sec. 14. II is indicated by 
the gapless rs-hu mode, while the gCFL" phase has both gs-bd (as in the 
gCFL phase) and rs-bu gapless modes. The definitions of all gapless phases 
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are summarized in the following table: 



Name 


Gapless mode(s) 
eik)^\k-kf\ 


Diquark condensate(s) 


g2SC 


ru-gd, gu-rd 


A3 


g2SC' 


rd-gu, gd-ru 


A3 


guSC 


rs-bu 


A2, A3 


gCFL 


gs-bd 


Ai, A2, A3 


gCFL' 


rs-bu 


Ai, A2, A3 


gCFL" 


gs-bd, rs-bu 


Ai, A2, A3 



4.4. Lepton fraction Yl^ 

Our numerical results for the lepton fraction Y^^ are shown in Fig. 1101 The 
thick and thin lines correspond to two different fixed values of the lepton- 
number chemical potential, /i^^ — 200 MeV and = 400 MeV, respec- 
tively. For a fixed value of /iL^, we find that the lepton fraction changes only 
slightly with temperature. This is concluded from the comparison of the re- 
sults at T = (solid hues), T = 20 MeV (dashed lines), and T = 40 MeV 
(dotted hues) in Fig. Uni 

As is easy to check, at T = 40 MeV, i.e., when Cooper pairing is not so 
strong, the /i dependence of Y^^ does not differ very much from the predic- 
tion in the simple two-flavor model in Sec. 14.11 By saying this, of course, 
one should not undermine the fact that the lepton fraction in Fig. 1101 has a 
visible structure in the dependence on at T = and T = 20 MeV. This 
indicates that quark Cooper pairing plays a nontrivial role in determining 
the value of Y^e- 

Our numerical study shows that it is hard to achieve values of the lepton 
fraction more than about 0.05 in the CFL phase. Gapless versions of the 
CFL phases, on the other hand, could accommodate a lepton fraction up to 
about 0.2 or so, provided the quark and lepton-number chemical potentials 
are sufficiently large. 

From Fig. 1101 we can also see that the value of the lepton fraction 
iie ~ 0-4: i-e., the value expected at the center of the protoneutron star 
right after its creation, requires the lepton-number chemical potential fi^^ 
in the range somewhere between 200 MeV and 400 MeV, or slightly higher. 
The larger the quark chemical potential ^ the larger /i^^ is needed. Then, 
in a realistic construction of a star, this is likely to result in a noticeable 
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Fig. 10. Dependence of the electron family lepton fraction Yi^^ for /i^^ = 200 MeV 
(thick lines) and /i^,^ = 400 MeV (thin lines) on the quark chemical potential at a fixed 
temperature, T = MeV (solid lines), T = 20 MeV (dashed lines), and T = 40 MeV 
(dotted lines). 

gradient of the lepton-number chemical potential at the initial time. This 
gradient may play an important role in the subsequent deleptonization due 
to neutrino diffusion through dense matter. 

4.5. T jjbL^ phase diagram 

Now let us discuss the phase diagram in the plane of temperature and 
lepton-number chemical potential, keeping the quark chemical potential 
fixed. Two such slices of the phase diagram are presented in Fig. ^2 The 
upper panel corresponds to a moderate value of the quark chemical poten- 
tial, — 400 MeV. This could be loosely termed as the "outer core" phase 
diagram. The lower panel in Fig. II II corresponds to ^ = 500 MeV, and we 
could associate it with the "inner core" case. Note, however, that the terms 
"inner core" and "outer core" should not be interpreted literally here. The 
central densities of (proto-)neutron stars are subject to large theoretical 
uncertainties and, thus, are not known very well. In the model at hand, the 
case /i — 400 MeV ( "outer core" ) corresponds to a range of densities around 
4/9o, while the case = 500 MeV ("inner core") corresponds to a range of 
densities around lOpo- These values are of the same order of magnitude 
that one typically obtains in models (see, e.g., Ref. 1^ . 

At first sight, the two diagrams in Fig. look so different that no obvi- 
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Fig. 11. The phase diagrams of neutral quark matter in the plane of temperature and 
lepton-number chemical potential at two fixed values of quark chemical potential: /i = 
400 MeV (upper panel) and fi = 500 MeV (lower panel). The triangle in the lower panel 
denotes the transition point from the CFL phase to the gCFL' phase at T = 0. 



ous connection between them could be made. It is natural to ask, therefore, 
how such a dramatic change could happen with increasing the value of the 
quark chemical potential from /i — 400 MeV to /z = 500 MeV. In order to 
understand this, it is useful to place the corresponding slices of the phase 
diagram in the three-dimensional diagram in Fig. 

The /i = 500 MeV diagram corresponds to the right-hand-side surface 
of the bounding box in Fig. This contains almost all complicated phases 
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with strange-nonstrange cross-flavor pairing. The /i — 400 MeV diagram, 
on the other hand, is obtained by cutting the three-dimensional diagram 
with a plane parallel to the bounding surface, but going through the middle 
of the diagram. This part of the diagram is dominated by the 2SC and NQ 
phases. Keeping in mind the general structure of the three-dimensional 
phase diagram, it is also not difficult to understand how the two diagrams 
in Fig. 111! transform into each other. 

Several comments are in order regarding the zero-temperature phase 
transition from the CFL to gCFL' phase, shown by a small black triangle 
in the phase diagram at /i = 500 MeV, see the lower panel in Fig. ^] The 
appearance of this transition is in agreement with the analytical result in 
Sec. 14.11 Moreover, the critical value of the lepton- number chemical poten- 
tial also turns out to be very close to the estimate in Eq. (|25() . Indeed, by 
taking into account that Mg « 214 MeV and A2 ~ 76 MeV, we obtain 
n'""^ = A2 — Mg/(2/i) « 30 MeV which agrees well with the numerical 
value. 

Before concluding this subsection, we should mention that a schematic 
version of the phase diagram in T-fj,L^ plane was earlier presented in Ref.|21 
see the left panel in Fig. 4 there. In Ref. |31 the value of the quark chemi- 
cal potential was /i = 460 MeV, and therefore a direct comparison is not 
straightforward. One can see, however, that the diagram of Ref. Elfits nat- 
urally into the three-dimensional diagram in Fig. Also, the diagram of 
Ref. i3 is topologically close to the ^ = 500 MeV phase diagram shown in 
the lower panel of Fig. ^2 The quantitative difference is not surprising: the 
region of the (g)CFL phase is considerably larger at /i = 500 MeV than at 
^ = 460 MeV. 

5. Conclusions 

Here we discussed the phase diagram of neutral three-flavor quark matter in 
the space of three parameters: temperature T, quark chemical potential /i, 
and lepton- number chemical potential ■ The analysis is performed in the 
mean-field approximation in the phenomenologically motivated NJL model 
introduced in Ref. |^ Constituent quark masses are treated as dynami- 
cally generated quantities. The overall structure of the three-dimensional 
phase diagram is summarized in Fig. and further detailed in several two- 
dimensional slices, see Figs. ITl 171 and ITTI 

By making use of simple model-independent arguments, as well as de- 
tailed numerical calculations in the framework of an NJL-type model, we 
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find that neutrino trapping helps Cooper pairing in the 2SC phase and 
suppresses the CFL phase. In essence, this is the consequence of satisfying 
the electric neutrality constraint in the quark system. In two-flavor quark 
matter, the (positive) lepton-number chemical potential helps to pro- 
vide extra electrons without inducing a large mismatch between the Fermi 
momenta of up and down quarks. With reducing the mismatch, of course, 
Cooper pairing gets stronger. This is in sharp contrast to the situation in 
the CFL phase of quark matter, which is neutral in the absence of elec- 
trons. Additional electrons due to large fi^^ can only put extra stress on 
the system. 

In application to protoneutron stars, the findings presented here suggest 
that the CFL phase is very unlikely to appear during the early stage of the 
stellar evolution before the deleptonization is completed. If color supercon- 
ductivity occurs there, the 2SC phase is the best candidate for the ground 
state. In view of this, it might be quite natural to suggest that matter in- 
side protoneutron stars contains little or no strangeness (just as the cores 
of the progenitor stars) during the early times of their evolution. In this 
connection, it is appropriate to recall that neutrino trapping also suppresses 
the appearance of strangeness in the form of hyperonic matter and kaon 

on 

condensation. The situation in quark matter, therefore, is a special case 
of a generic property. 

After the deleptonization occurs, it is possible that the ground state of 
matter at high density in the central region of the star is the CFL phase. 
This phase contains a large number of strange quarks. Therefore, an abun- 
dant production of strangeness should happen right after the deleptoniza- 
tion in the protoneutron star. If realized in nature, in principle this scenario 
may have observational signatures. 
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